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We design linear optics multiqubit quantum logic gates. We assume the traditional encoding 
of a qubit onto state of a single photon in two modes (e.g. spatial or polarization). We suggest 
schemes allowing direct probabilistic realization of the fundamental Toffoli and Fredkin gates without 
resorting to a sequence of single- and two-qubit gates. This yields more compact schemes and 
potentially reduces the number of ancilla photons. The proposed setups involve passive linear 
optics, sources of auxiliary single photons or maximally entangled pairs of photons, and single- 
photon detectors. In particular, we propose an interferometric implementation of the Toffoli gate 
in the coincidence basis, which does not require any ancilla photons and is experimentally feasible 
with current technology. 
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I. INTRODUCTION 

Quantum information theory Q exploits the laws of 
quantum mechanics to devise novel methods of informa- 
tion processing and transmission that would be impos- 
sible or very hard to achieve classically. During recent 
years various protocols for quantum information process- 
ing were successfully demonstrated experimentally with 
several different physical systems. Particular attention 
has been paid to optical realizations where the quantum 
bits are encoded onto states of single photons. Photons 
are ideal carriers of quantum information because they 
can be distributed over long distances in low-loss optical 
fibers or in free space. While perfect for quantum com- 
munication purposes, photons seemed to be less suitable 
for quantum computing because the lack of sufficiently 
strong optical nonlinearities seemed to prevent the im- 
plementation of quantum gates between photons. 

The situation changed radically in 2001 when Knill, 
Laflamme, and Milburn (KLM) published their landmark 
paper in which they showed that a scalable universal 
quantum computation is possible with only single photon 
sources, passive linear optical interferometers and single 
photon detectors 0|- The key insight of KLM is that 
the nonlinearity (such as a Kerr effect) can be simulated 
on a single-photon level using the above listed resources, 
conditioning on particular measurement outcomes of the 
detectors and applying appropriate feedback. The re- 
sulting linear optics quantum gates Q are generally only 
probabilistic but the probability of success could be in 
principle made arbitrary close to unity by exploiting off- 
line generated multi-photon entangled states and quan- 
tum teleportation ^ ||6i . 

The KLM paper stimulated a number of further works 
suggesting alternative and improved constructions of the 
basic quantum C-NOT gate B II H IS IH El whose 
experimental demonstrations by several groups followed 
[IlllllllllllllllllliliSlill- However, despite 
these promising successes, extending this approach to 
more complex schemes involving higher number of pho- 
tons currently appears to be a formidable experimental 



task because the overhead in resources (in particular the 
number of ancilla photons) required by the original KLM 
scheme is very high. 

It is possible to combine the ideas of one-way quantum 
computation (2% . and linear optics quantum comput- 
ing to significantly reduce the resources required for the 
computation. The techniques introduced by KLM could 
be used to generate a multiphoton cluster state which 
then serves as a resource for quantum computing which 
proceeds by performing certain carefully chosen measure- 
ments on each photon from the cluster 24, 25, 26, 27|. 
First proof-of-principle experimental demonstration of 
one-way quantum computation with four-photon cluster 
state has been reported recently [28]. 

In this paper, we wish to address a different aspect of 
the quantum computing with linear optics. Namely, we 
will be interested in the implementations of the funda- 
mental Toffoli and Fredkin gate which play an important 
role both in classical (reversible) computing and in quan- 
tum computing and information processing P|. In the 
universal quantum computer the multi-qubit gates are 
usually assumed to be implemented as a sequence of sin- 
gle and two-qubit gates. However, this strategy may not 
be optimal in the context of linear optics quantum com- 
puting, where schemes tailored specifically for multiqubit 
gates may require less ancilla photons or achieve higher 
probability of success than implementations relying on a 
sequence of the single and two-qubit gates. 

The rest of the paper is organized as follows. In Section 
II we will present a scheme for the iV-qubit generalized 
Toffoli gate, which flips the state of the A^th qubit in the 
computational basis if all the — 1 control qubits are 
in state |1), where tilde indicates the logical qubit states 
throughout the paper to distinguish them from the Fock 
states \n). We will first design gate operating in the 
so-called coincidence basis which has the advantage 
that this scheme does not require any ancilla photons. To 
make this gate non-destructive it is necessary to perform 
quantum non-demolition measurement of number of pho- 
tons at the output of the gate, which could be done with 
linear optics, ancilla photons and photon-number resolv- 
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FIG. 1: Setup for the linear-optics A^-qubit controlled phase 
gate. Tj label beam splitters with transmittance Tj, M in- 
dicate mirrors, PS is a phase shift by amount i/;. Each grey 
ellipse represents a single qubit carried by a single photon 
propagating in two different spatial modes. 

ing detectors |30|. Section III is devoted to the three- 
qubit Fredkin gate which is a controlled SWAP gate, the 
states of the two target qubits are swapped if the control 
qubit is in state Our scheme requires only six ancilla 
photons. Finally, Section IV contains brief conclusions 
and summary of the main results. 



II. QUANTUM TOFFOLI GATE WITH LINEAR 
OPTICS 

In this section we will present and analyze the scheme 
which realizes quantum iV-qubit Toffoli gate in the coin- 
cidence basis. More precisely, the scheme conditionally 
applies the A^-qubit controlled phase gate to the input 
qubits, whereby the phase changes by tt if all qubits are 
in the state |1) and does not change otherwise, 

C/SpIji, J2, . . . , Jiv) = e^"n:.i:"^-|j-i, . . . , j^), (1) 

where jk € {0,1} and k = 1,...,N. Note that the 
N-qubit controUed-phase (C-phase) gate and the Toffoli 
gate are equivalent up to single-qubit Hadamard trans- 
formations H on the A''th qubit (the target), Ut ~ 
HnUcpHn. 



A. Quantum optical C-phase gate 

The proposed optical setup is schematically sketched in 
Fig. 1. The qubits are encoded into states of single pho- 
tons. We assume the dual-rail encoding where the two 
logical levels \0)j and correspond to two paths Jl and 
jfl: taken by a photon, \0)j = \01)j^j^ and \l)j = 110)^^^^. 
Note that other encodings such as polarization or time- 
bin are also possible and can be mutually converted into 



each other by means of polarizing beam splitters and un- 
balanced interferometers. 

The operation of the C-phase gate requires that the 
phase of the iV-photon state changes by tt if and only if 
all photons are in the L modes. In the proposed scheme 
this is achieved by the A^-photon interference on an array 
of unbalanced beam splitters, see Fig. 1. First, each pho- 
ton propagating in the L mode is split into two modes 
on a beam splitter with intensity transmittance Ti . Then 
pairs of beams originating from modes jl and (j -I- 1)l 
interfere on an array of N beam splitters with transmit- 
tance T2. Each mode passes through a beam splitter 
with transmittance T3 which acts as a filter balancing the 
amplitudes of the modes Jl and jn after the application 
of the gate. 

The gate operates in the coincidence basis, i.e. it suc- 
ceeds if a single photon is detected in each pair of the 
output modes Jl and ju. Assume first that at least one 
photon is in mode fc^j. It is easy to see that in this case 
the only way how the photons in the ji modes can reach 
the appropriate output ports of the beam splitters is that 
each photon is transmitted through both beam splitters 
and the total probability amplitude of this to happen 
reads a„ = (tii2)^~"i3, where n > 1 is the number of 
photons in kfj modes and tj is the amplitude transmit- 
tance, Tj = t'j. Since the gate should be unitary, a„ must 
not depend on n which can be achieved by choosing 

T3 = T1T2. (2) 

The situation changes when all photons are initially 
in L modes, i.e. the input state reads |1, 1, 1, ... , 1). In 
this case, there are two ways how the photons can reach 
the A^ output ports. One option is that all photons are 
transmitted through all beam splitters. The second op- 
tion is is that all photons are reflected from all beam 
splitters. Provided that these two alternatives are indis- 
tinguishable, i.e. there is a good spatiotemporal overlap 
of the photonic wavepackets on the beam splitters, they 
interfere, and the resulting amplitude reads, 

ao = «-rfrf. (3) 

Note the minus sign which arises due to the tt phase 
shift 1/) = TT in one path of the reflected photon, see Fig. 
1. Note also that in the figure the path of the refiected 
photon in mode looks much longer than all other 
paths. In the actual implementation the geometry of the 
setup should be such that all paths would be carefully 
balanced resulting in a good overlap of the photons and 
high- visibility interference. 

The gate operates as desired if ag = — a„>o — — t^i^. 
Expressed in terms of the intensity transmittances this 
condition translates into 

AT^T^ = (1 - ri)^(l - Ts)^, (4) 

where we used that Vj ~ y/1 — Tj. The formula Q de- 
scribes a single-parametric class of the N-qubit optical 
controUed-phase gates working in the coincidence basis. 
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The probability of success is given by Psucc = = 
j'NrpN g^j^j expressing T2 in terms of Ti from Eq. Q), 



T2 



l-Ti 



we obtain 



1 - Ti + ' 
Ti(l-Ti) 



1 - Ti + 



N 



(5) 



(6) 



The optimal Ti maximizing Pgucc can be easily deter- 
mined by solving dPgucc/dTi — 0, which yields 



Ti.o- 



1 



1 + 21/w 



(7) 



On inserting this value into Eq. @ we find that r2_opt = 
T'i,opt hence it is optimal to use a scheme where the trans- 
mittances Ti and T2 are the same. The optimal proba- 
bility then reads, 



succ.opt 



1 



1 + 



2N 



(8) 



In particular, for = 3 we find P, 



N=3 
succ,opt 



0.75%. 



B. Generalized C- phase gate 

The transformation Q can be extended such that an 
arbitrary phase shift (j) is introduced when all qubits are 
in logical state The generalized controUed-phase gate 
thus acts as follows, 

U^p\ji,j2, ■■■,Jn)^ e^^nlLi^l^-i, . . (9) 

We shall show that also this operation can be condition- 
ally implemented with the scheme shown in Fig. 1 pro- 
vided that the phase shift tp in one arm of the multipho- 
ton interferometer and the transmittances Ti and T2 are 
properly chosen. Repeating the derivation outlined in 
the preceding subsection we find that the condition that 
has to be satisfied reads 



(10) 



Upon splitting this formula into the real and imaginary 
parts and solving for tp we obtain 



tan?/' = — - 



1 



tan ^ 



(11) 



Similarly we also arrive at a generalization of the formula 
0, 



4 sin^ I = (1 - Ti)^ (1 - T^) 



N 



(12) 



The probability of success is maximized by choosing 



l-f |2sin Ill/A' ^ 



(13) 
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FIG. 2: Equivalence between two-qubit controlled unitary 
and controlled phase gates. 



and we have Pc 



succ.opt 



(0) = 1/[1 + |2sin((/)/2)|i/^]2^. 
Note that the probability of success depends on the re- 
quired conditional phase shift (jj and for a fixed N it 
achieves its minimum for (jj — tt, i.e. when we attempt to 
implement the N-qubit Toffoli gate. 

With the two-qubit generalized C-phase gate at hand 
we can implement in the coincidence basis an arbitrary 
two-qubit controlled-J7 gate, where a unitary operation 
U is applied to the target qubit iff the control qubit is 
in state |l)c- The equivalent scheme involving C-phase 
gate is shown in Fig. 2(b). Note first that in the basis of 
eigenstates |uj)T of C/, the controlled unitary gate boils 
down to conditional phase shifts, 

\0)c\uj)t ^ \0)c\uj)t, \i)c\uj)T ^ e'*^-|i)c|u,)T. 

(14) 

The unitary V maps the eigenstates of U onto the com- 
putational basis states, V\uj)T = \j)T, j = 0, 1. Next 
a C-phase gate with A0 = 0i — 0o follows. Finally, 
the inverse operation is applied to the target while 
the control qubit is subject to a phase shift operation 
|0)c \0)c, |i)c e*'^°|i)c. It is easy to see that the 
net result of this sequence of gates is the controlled- C/ 
operation 114|1 . 



C. Heralded controlled-phase gate 

The advantage of working in the coincidence basis is 
that no extra ancillary photons are required. The scheme 
is thus very economical in resources and for instance the 
demonstration of the three-qubit Toffoli gate would re- 
quire detection of three-photon coincidences which is well 
within the scope of present technology. 

However, this approach also suffers from a significant 
disadvantage, since we do not know whether the gate 
succeeded until we detect the photons. It is thus not 
possible to directly employ this gate as a part of a more 
complex quantum information processing network. Nev- 
ertheless, it is possible to remove this drawback by per- 
forming quantum non-demolition measurements of the 
number of photons at the outputs of the gate. If this 
measurement verifies that a single photon is present in 
each pair of modes and jn then we know that the 
gate was applied successfully while the non-demolition 
character of this measurement guarantees that the out- 
put photons emerging from the gate are preserved and 
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not destroyed by the verification. 

A simple way of performing the non-demoUtion mea- 
surement of a number of photons in two modes is to cm- 
ploy an auxiliary pair of photons in a maximally entan- 
gled state and attempt to teleport the single photon in 
modes Jl and Jr [s^I • The single-photon detectors used 
for the partial Bell measurement which lies at the heart of 
the teleportation j5i| must be able to resolve the number 
of photons. Detection of exactly two photons in the Bell 
analysis confirms that a single photon has been success- 
fully teleported. Observation of any other total number 
of photons indicates a failure of the gate. This method 
requires N auxiliary maximally entangled photon pairs 
in total and the probability of successful non-demolition 
measurement given that the gate was applied successfully 
scales as 1/2^ because the optimal partial Bell measure- 
ment with linear optics can distinguish only two out of 
four Bell states. 

An alternative scheme for partial probabilistic non- 
demolition photon number measurement on a pair of 
modes has been proposed in 's^ (see also discussion in 
Rcf. 31] ). The advantage of this latter scheme is that it 
does not rely on maximally entangled photon pairs and 
instead requires single photons in product state, which 
may be easier to generate. The measurement requires 
two ancilla photons and two photodetectors which can 
distinguish the number of photons in a mode. A co- 
incidence detection of a single photon by each detector 
indicates that at least a single photon has been present 
in the input pair of modes and if exactly a single photon 
was at the input then its state was not disturbed by the 
measurement. If this partial measurement is carried out 
on each pair of modes , ja and if all N measurements 
indicate that there was at least a single photon in each 
pair of modes then since there were altogether N pho- 
tons at the input of the gate we can conclude that the 
C-phase gate was applied successfully. 



III. FREDKIN GATE 

Our scheme for linear optics Fredkin gate is inspired 
by the quantum optical Fredkin gate originally proposed 
by Milburn (s^l- The Fredkin gate is a controlled SWAP 
operating on the Hilbert space of three qubits, the states 
of qubits A and B are exchanged if the control qubit C is 
in state |1} and nothing happens if it is in state |0}. Let 
us assume that the qubits are encoded onto polarization 
states of single photons. The controlled SWAP operation 
can be converted to the controlled phase shift with the 
use of a balanced Mach-Zehnder interferometer, see Fig. 
3. Depending on the state of the control photon C, the 
phase shift in the left arm of the interferometer should be 
either or tt, the latter results in the effective swap of the 
photons A and B at the output of the interferometer. In 
Milburn's scheme, the controlled phase shift is achieved 
by medium with cross Kerr nonlinearity. 

In the spirit of linear optics quantum computing, we 




FIG. 3: Quantum optical Fredkin gate based on a bal- 
anced Mach-Zehnder interferometer. BS denote balanced 
beam splitters and M indicate mirrors. If the control pho- 
ton in mode C is vertically polarized then the boxes CPSV 
and CPSH apply conditional phase shift tt to the vertically 
(CPSV) or horizontally (CPSH) polarized modes in the left 
arm of the interferometer. 

suggest to replace the Kerr medium with a linear interfer- 
ometric scheme, and employ ancilla photons and postse- 
lection conditioned on single photon detection to simulate 
the required cross Kerr interaction ■ Since we assume 
polarization encoding, the conditional phase shift has to 
be applied to both vertically and horizontally polarized 
modes in the left arm of the Mach-Zehnder interferome- 
ter in Fig. 3. In what follows we will describe a scheme 
which provides this conditional phase shift for a single 
mode, and the linear optics Fredkin gate then involves 
two such basic blocks acting in series on vertically and 
horizontally polarized mode, see Fig. 3. 

The basic block is depicted in detail in Fig. 4. The 
scheme requires three ancilla photons: a maximally en- 
tangled pair of photons in a state ( I + 1 HH) ) emit- 
ted by source (EPR) and an additional single photon in 
mode 2. The proposed setup consists of two main parts. 
The first part is the quantum parity check between 
the control photon in mode C and one photon from the 
auxiliary EPR beam. The check is based on a coupling 
of these photons on a polarizing beam splitter PBS fol- 
lowed by a detection of one of the outputs in the basis 
^(|y) ± \H)). The detectors should be able to resolve 
the number of photons in the beam and the parity check 
is successful if a single photon is detected by one de- 
tector and no photon is detected by the other detector, 
which happens with probability 1/2. The parity check 
effectively copies in the computational basis the state of 
the control photon in spatial mode C onto the auxiliary 
photon in spatial modes 3 and 4 and we can write, 

a\V)c + l3\H)c ^ a\V)c\0)3\l)4 + m)c\l)3\0)4. (15) 

The parity check allows us to control the phase shift of 
mode 1 indirectly by the auxiliary photon while preserv- 
ing the original control photon. A similar trick has been 
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FIG. 4: Simulation of cross-Kerr interaction with linear op- 
tics. The setup involves polarizing beam splitters (PBS), 
source of an auxiliary entangled pair of photons (EPR) , source 
of a single photon (|1)), linear multiport interferometer con- 
sisting of beam splitters and phase shifters and single-photon 
detectors (PD). 



used in the recent experimental implementations of quan- 
tum C-NOT gate .14, 17, 18]. 

The second part of the scheme in Fig. 4 consists of 
a linear interferometer where the photons in the mode 1 
are combined with the ancilla photons in modes 2, 3 and 
4. Note that the interferometer has also three other aux- 



iliary input ports in vacuum state. All output modes of 
the interferometer except for mode 1 are monitored with 
photon number resolving detectors and the conditional 
phase shift is successfully applied if a single photon is 
detected in modes 2 and 3 and no photons are observed 
in the other modes. 

The purpose of the interferometer is to conditionally 
induce a phase shift tt in mode 1 provided that there 
is a photon in the input mode 4 and induce no shift if 
the photon is in mode 3. Since there can be no more 
than 2 photons in the mode 1 (the two photons whose 
polarization states should be conditionally swapped), it 
suffices to achieve the correct conditional phase shift in 
the subspace of Fock states |0)i, |l)i and |2)i. 

Mathematically, the interferometer is described by a 
unitary matrix U , which governs the transformation be- 
tween input and output modes. The input creation op- 
erators ajjj J can be expressed as linear superpositions of 
the output creation operators a^^j j. according to 



fc=i 



(16) 



where Ujk are the elements of U. Since we condition on 
observing no photons in modes 4-7, in our subsequent 
calculations we will explicitly need only the coefficients 
Ujk with j = 1, 2, 3, 4 and fc = 1, 2, 3. 

Due to the linearity we can treat separately the cases 
when the control photon is in state \H) and \ V). Assume 
first that it is in state \H). The two auxiliary photons 
are then in input modes 2 and 3 and conditionally on 
detecting a single photon in the output modes 2 and 3 
and no photon is all other modes (except for mode 1 
which is not measured upon), we obtain the following 
transformation. 



|if)c|l)2|l)3|0)4(ao|0)i + ai|l)i + a2|2)i) ^ |i/)c(a;oao|0)i + 2;iai|l)i + X2a2\2)i) (17) 
where the coefficients xj can be expressed as follows, 

Xq = U22"33 + U23U32, 

Xi = Ull(u22"33 + U23W32) + M12(U23W31 + W21U33) + U13(U22U31 + M21M32), 

X2 = U?i(w33"22 + M32W23) + 2M12U13U21U3I + 2W11W12(U23W31 + U21U33) + 2U11M13(U22'"31 + U21"32)- (18) 

If the control photon is in state \V), then the conditional transformation reads 

|^^)c|l)2|0)3|l)4(ao|0}i + ai|l}i -I- a2|2)i) |T^)c(yoao|0)i + yiai|l)i + y2a2|2}i), (19) 



where the coefficients yj can be expressed in the same 
way as Xj, only the matrix elements u^k iu Eq- l|18() 
must be replaced with U4k. 

We want to implement a conditional 7r-phase shift in 



mode 1. This will be achieved if 

x,^q, y,=q{-iy, (20) 
where q < 1 is some shrinking factor arising due to the 
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probabilistic nature of the gate. Low q reduces the prob- 
abihty of success of the gate which scales as P oc \q\'^ 
but does not alter its operation. The maximum q that 
can be chosen is determined by the constraint that Ujk 
(1 < j < 4, 1 < A: < 3) must form a submatrix of a 
unitary matrix. As shown in the Appendix, it is possi- 
ble to efficiently numerically determine whether a given 
set of Ujk may form a submatrix of U so that also the 
maximum q can be determined numerically. 

By solving the system of nonlinear equations H2U|) we 
can find the matrix elements Ujk specifying the inter- 
ferometer which implements the conditional phase shift. 
Note that the system is underdetermined hence there ex- 
ist infinitely many interferometers satisfying (|20|1 . To see 
this, it is convenient to rewrite these equations in a ma- 
trix form, 




Mua 




(21) 



where 1x3 = (uai, "32, ""33)^ and U4 = (u4i, M42, ^43)^ are 
column vectors and M is a matrix whose elements can 
be expressed in terms of uik and U2k- Thus when looking 
for the solution to Eqs. H21|l we can choose arbitrary uik 
and U2k and provided that det M 7^ we can for a given 
q calculate 1x3 and 1x4 by solving the system of linear 
equations (|^ . 

Let us now present a particular example of an analyt- 
ical solution. Choosing 



un = U12 = U13, 



"21 



'U22 = "23, 



(22) 



we can express the matrix elements u^k as follows, 



U31 
U32 
U33 



2UiiU22 



2uiiU22 ' 



(2uii ~ 1). 



2uii?i22 

Similar formulas hold also for u^k and we have, 



(23) 



U41 
U42 
U43 



2u22Un 
g 

2wiiU22 ' 



2wilU22 



(2uii + l). 



(24) 



The maximum \q\ achievable within the above given 
analytical solution was determined numerically and we 
found that it is optimum to choose wn.opt = 0.494 and 
M22,opt — 0.416 yielding gopt — 0.0638. Since the Fredkin 
gate in Fig. 3 includes two conditional phase shift gates, 
the total probability of success of the gate is given by 
■Psucc = j\q\'^ where the factor j appears due to the two 
quantum parity checks. On inserting the (/opt into this 



formula we obtain P„ 



4.2 X 10 ^, which is rather 



small. However, it should be stressed that this is not 
the maximum probability of success that could be at- 
tained with our scheme. It is possible to improve the 
success rate by several orders of magnitude by perform- 
ing numerical optimization over all relevant parameters 
uik and U2k- We have carried out a thorough numerical 
search and the maximum Psucc that we obtained in this 
way reads Psucc,max = 4.1 x 10~^. 

In is instructive to compare this value with the prob- 
ability of success that could be achieved if one would 
attempt to implement the Fredkin gate as a sequence of 
two-qubit unitaries. It was shown by Smolin and DiVin- 
cenzo that five two-qubit quantum gates suffice to imple- 
ment the Fredkin gate |3J|. Making the very optimistic 
assumption that using two ancilla photons per gate each 
of these gates can be irnplemented with probability 1/4 
similarly as the C-NOT '^^ arrive at a total prob- 

ability Pl^^^ = w 9.8 X 10"''. Thus our scheme for 
Fredkin gate could potentially attain a higher probabil- 
ity of success while being more economical in resources 
because it requires only 6 ancilla photons instead of 10 
photons. 



IV. CONCLUSIONS 

We have devised schemes for linear optics quantum 
Toffoli and Fredkin gates. In the spirit of linear optics 
quantum computing the gates do not require nonlinear 
interaction and instead rely on multiphoton interference, 
ancilla photons and postselection conditioned on single- 
photon detection. The key feature of the proposed setups 
is that they are directly tailored for the implementation 
of the multiqubit Toffoli or Fredkin gate. This should 
be contrasted with the common approaches where the 
multiqubit gates are decomposed into a sequence of two- 
and single-qubit gates. 

Given the current state of the technology, our direct 
approach to multiqubit gates may be much more efficient 
than implementations relying on a sequence of two-qubit 
gates. In particular, the experimental demonstration of 
the three-qubit quantum Toffoli gate in the coincidence 
basis would require only three photons and an observa- 
tion of a three-photon coincidences, which is well within 
the reach of current technology. 

Despite their advantages, the present schemes still suf- 
fer from some weaknesses. The probability of success 
of the N-qubit Toffoli gate exponentially decreases with 
growing N and also the probability of success of the Fred- 
kin gate, Psucc ~ 4.1 X lO^'^, is not very high. Another 
drawback lies in the fact that setups for both Toffoli and 
Fredkin gate require interferometric stability, which is 
hard to achieve and maintain. In contrast, recent exper- 
imental demonstrations of the quantum linear-optical C- 
NOT gate relied solely on Hong-Ou-Mandel interference 
effect |l9l I20I l2ll |. which is much more robust against 
small length fluctuations. 

It remains an interesting open question whether a 
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scheme avoiding problems with interferometric stabihty 
could be devised also for the Toffoli and Fredkin gates. 
Another important open issue is what is the maximum 
achievable success rate for these gates either without an- 
cillas (i.e. operating in the coincidence basis) or with a 
given fixed amount of auxiliary photons. We hope that 
the present paper will stimulate further theoretical as 
well as experimental investigations along these lines po- 
tentially resulting in an important step towards linear 
optics quantum computing. 
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APPENDIX: TESTING THE UNITARITY 



ative and determine it from the normalization condition 



— 



\ 



J+2 
k=l 



If this expression yields purely imaginary then Ujk 

could not form a part of a unitary matrix and we termi- 
nate the test. 

Step (ii): We must guarantee that the rows of U are 
mutually orthogonal. This can be achieved by calculating 
Uk.j+3 from the orthogonality condition 



E u*iUki = 0, 



For all k satisfying j < k < A we thus have 



Uk,j+3 



J+2 



Here we show how to test whether a 4 x 3 matrix ujk , 
1 < J < 4, 1 < A: < 3, could be a submatrix of a larger 
7x7 unitary matrix U. The procedure consists of two 
steps which have to be repeated for each j e {1,2,3,4}. 
We start from j — 1. 

Step (i): We set ujk = for fc > j + 3. We can assume 
without loss of generality that itj.j+3 is real and nonneg- 



Finally, we increase j by 1. 

The steps (i) and (ii) have to be repeated until j — 
4 is reached. If all four iterations succeed then at the 
end we obtain a 4 x 7 isometry matrix that could be 
easily completed to form a unitary matrix. Otherwise we 
know from step (i) that the matrix Ujk could not form a 
submatrix of a unitary matrix. 
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